Optimization of Short Coherent Control Pulses 



00 

o 
o 

(N 



^ ■ 

a ■ 



> 
oo 

00 

p 

O 
O 



X 



S. Pasini/ T. Fischer/ P. Karbach/ and G. S. Uhrig^ 

' Lehrstuhl fur Theoretische Physik I, Technische Universitdt Dortmund, 
Otto-Hahn Strafie 4, 44 Dortmund, Germany 
(Dated: March 11, 2008) 

The coherent control of small quantum system is considered. For a two-level system coupled to 
an arbitrary bath we consider a pulse of finite duration. We derive the leading and the next-leading 
order corrections to the evolution operator due to the non-commutation of the pulse and the bath 
Hamiltonian. The conditions are computed that make the leading corrections vanish. The pulse 
shapes optimized in this way are given for n and ^ pulses. 
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I. INTRODUCTION 

For a long time, the coherent control of quantum sys- 
tems has been an important issue. For instance, nuclear 
magnetic resonance (NMR) is a widely used technique 
for the analysis of chemical compounds and for imag- 
ing techniques. Another broad field of application is the 
coherent control of the state of quantum dots, both in 
the charge and in the spin degrees of freedom. In recent 
years, however, coherent control has gained an even more 
vivid interest because it is the indispensable prerequisite 
in quantum information processing. 

In its most basic form, the quantum state of a single 
quantum bit (qubit) is to be changed in a given way. This 
is called a single qubit gate. A qubit is a two-level system 
which can be conveniently seen as S = 1/2 system. So we 
use this spin language to describe the qubit gates. The 
two most common gates are the tt pulse, which flips the 
spin from up to down and vice- versa, and the 7r/2 pulse 
or Hadamard gate, which rotates the spin by 90° away 
from the Sz direction. 

The TT pulse isparticularly interesting for dynamic de- 
couphng [I!, Hi, [111, [1,0, 111- In a nutshell, this technique 
aims at decoupling the single qubit from its environment 
as well as possible by switching the qubit state by single 
qubit gates. Mostly, the tt pulse is considered. The idea 
comes from the spin-echo technique in NMR where a tt 
pulse is used to refocus the precessing magnetization. In 
NMR, the use of sophisticated protocols of tt and 7r/2 
pulses is common to suppress unwanted couplings be- 
tween nuclear spins [§, [l3| • 

Experiments using tt pulses have obtained many en- 
couraging results in prolonging the coherence time of a 
qubit [ill, [13, [H, [iJI- More complex pulse sequences 
based on 7r/2 pulses have also proved useful by deminish- 
ing the rate of decoherence in quantum registers consid- 
erably [H. 

In the theoretical treatments the generic ingredient 
to dynamic decoupling is an idealized pulse, mostly a 
TT pulse, which is considered to be instantaneous. This 
means that its amplitude is infinite in the sense of a 
(5-function. This assumption is convenient because the 
Hamiltonian Hq of the pulse and the Hamiltonian H 



of the system at rest do not commute in general. The 
Hamiltonian H comprises the Hamiltonian of the qubit 
itself, the Hamiltonian of the environment (henceforth 
called "bath") and the coupling between the qubit and 
the bath which is responsible for its decoherence. But 
during the pulse Hq dominates and H is negligible and 
for the rest of the time only H is present. Hence it is rel- 
atively straightforward to deal with the time evolution. 

Experimentally, however, this idealized situation is not 
realistic. The real pulses are always of finite amplitude 
and of finite length. So the question has to be addressed 
which effect a finite pulse duration Tp has. We want to 
elucidate this issue in the present work theoretically in a 
fully quantum mechanical framework. 

The issue of pulses of finite duration and finite ampli- 
tude has been considered in other papers. Viola and Knill 
generalized their previous approach of averaging over a 
symmetry group [l6| by instantaneous changes of the ef- 
fective Hamiltonian to continuous ch ang es requiring only 
finite, bounded control amplitudes [l7|- But the pres- 
ence of the coupling to the bath during the pulses has 
not been considered. 

Khodjasteh and Lidar considered the effect of finite 
pulse duration in a comprehensive comparison of various 
schemes for n pulse sequences [11]. They did not con- 
sider to make the detrimental effects of the finite pulse 
duration vanish on the level of the individual pulse, but 
they discussed possible cancellations on the higher hier- 
archical level of the pulse sequence. In the present work, 
we will choose another route and aim at reducing the 
decoherence due to finite pulse duration already on the 
level of the single pulse. As a side remark, we note that 
the statement in Ref. [l^ on optimized pulse sequences 
S] that they hold only for bosonic bath operators is very 
likely not to be correct. It has been conjectured very 
recently by Lee et al. [l^ that the optimized sequences 
apply to the most general phase decoherence model. 

Very recently, the effect of classical random telegraph 
noise during pulses of finite length has been discussed 
numerically [20|. It was shown that shaping the pulse 
amplitude in a particular way can be used to improve 
its performance. This means that one can optimize the 
pulse such that its effect is closer to the desired one, e.g. 



a TT pulse. 

Tuning pulses by shaping them is a possibility of opti- 
mization which has been discussed intensively in the vast 
literature on NMR. But the objectives were mostly dif- 
ferent from the consideration of the effects of a bath. For 
instance, Geen and Freeman aimed at frequency selectiv- 
ity i.e. the pulse should act only on spins on or close 
to resonance, but leave others unchanged. Closest to our 
scope are the investigations by Tycko [23 |. Cummins et 
al. [23I . [23 | , and Brown et al. [2^ where composite pulses 
are optimized to be robust against off-resonance effects 
and pulse length inaccuracies. Again, no dynamic but 
only static effects are considered. Interestingly, we will 
see that some of the pulse shapes fulfilling the conditions 
in first order derived here for fully dynamical baths coin- 
cide with pulse shapes previously derived to compensate 
static effects. 

In this work, we present a systematic expansion in the 
pulse duration Tp considering a qubit coupled to a fully 
quantum mechanical and dynamical bath. This means 
that the related energy scale LUp :— l/rp (h is set to unity 
for simplicity) is considered to be very large so that all 
other energy scales, namely the coupling A between qubit 
and bath and the internal energy scale of the bath ivb, 
are small relative to Up. We expand in X/tOp — TpX and 
LUii/ujp =^ TpLOf). The zeroth order is the instantaneous 
pulse. We compute the first and the second order and 
discuss the conditions for which they vanish. These con- 
ditions allow us to predict pulse shapes which approx- 
imate instantaneous pulses in spite of their finite pulse 
duration. 

We emphasize that it is not our primary aim to elimi- 
nate the coupling between qubit and bath. On the con- 
trary, it turns out that the expansion around the instan- 
taneous pulse keeps this coupling. But the coupling to 
the bath is disentangled from the actual pulse. They are 
expressed by separate, subsequent time evolution oper- 
ators. Given the thus optimized pulse shapes, the cou- 
pling to the bath can be compensated on the higher hi- 
erarchical level of an appropriate pulse sequence, i.e., by 
dynamic decoupling y, H, d, H, [Ti, i8j] . 

The paper is organized as follows. After introducing 
the model (Sect. II), we start describing the method we 
used to derive our general equations for the first and sec- 
ond order (Sect. III). Then we discuss the results and 
provide various examples for pulses that satisfy our con- 
ditions (Sect. IV). At last we discuss our results and com- 
pare them with other proposals for shaped pulses in the 
literature (Sect. V). 



II. MODEL 

We consider the time evolution operator from t = to 

t = Tp 

Up{Tp,0) = Tl^cxp(^-i J^' Ho{t)dt-iHTp^^ . (1) 
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Figure 1: (a) Plot of a generic pulse v{t). The x- and y-axis 
values are expressed in units of Tp and l/rp, respectively. The 
shaded region represents an example for 'p{t), that is the inte- 
grated pulse amplitude between the instant and a generic 
instant t (see text). 

(b) Sketch of the approximately equivalent idealized instan- 
taneous pulse. 

For simplicity the pulse is chosen to start at i — 0; T 
stands for the usual time-ordering. The control Hamil- 
tonian of the pulse shall be given by 

Ho^v{t)ay, (2) 

where v{t) stands for the pulse shape as function of time, 
see Fig. [TJ The time-independent Hamiltonian of the 
bath and the coupling to the qubit reads 

H = Hb + \Aa,, (3) 

with Hfj representing a generic bath and A its coupling 
operator to the qubit. The internal energy scale of Hh 
shall be denoted by ujb- Note that Hq and H do not 
commute so that the evolution operator ([1]) is a non- 
trivial expression. 

The model ([3]) does not include spin flips; hence it im- 
plies an infinite life time Ti. But the decoherence of a 
precessing spin in the xy plane is described in full gen- 
erality because we do not specify for which operator A 
stands and the bath dynamics described by Hi,. Such a 
model is experimentally very well justified as the effec- 
tive model in the limit of a large applied magnetic field 
which implies that other couplings between the quantum 
bit spin and the bath are averaged out, sec for instance 

Ref . 0, Ha. 

As an example only one may consider the spin-boson 
model 

H^Y. + '^^ E ^'(^« + (4) 

i i 

with A := maxiAj, A = J2ii^i/^)i^l + ^i) ^"^'^ ■= 
maxi LUi . 
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The primitive of the pulse amphtude is 



(5) 



where the integration starts from an instant Tg G [0,Tp]. 
The intended angle of rotation is given by twice the total 
area under the curve 9 :— 2^ — 2{(p{Tp) — ip{Oj)- This 
area has to be kept constant while investigating the limit 
Tp ^ so that the amplitude v{t) has to tend to infinity 
as v{t) oc 1/Tp. To simplify our notation later, we define 
also ip{t) 2ip{t). 

The aim is to replace the real pulse given by v{t) by 
an approximately equivalent instantaneous pulse acting 
at the instant r^, see Fig. [T] The equivalence shall hold 
up to linear order or quadratic order in Tp. We highlight 
that the time evolution before and after the equivalent in- 
stantaneous pulse shall be governed by H, i.e., including 
the decohering interaction between qubit and bath. This 
shall be accounted for on the higher hierarchical level of 
appropriate pulse sequences. 



III. METHOD AND GENERAL EQUATIONS 

A. General Equation 

Let Ts be an arbitrary instant inside the interval [0, Tp] 
at which we want to split the time evolution 



Up{Tp, 0) = Up {Tp, Ts) Up (Ts, 0) , 



(6) 



The guiding idea is shown in Fig.[T]3. We want to replace 
the full evolution by a first part during [0, r^] where only 
the Hamiltonian H is active. Then the instantaneous 
pulse shall follow. Finally, another intervall [rs,Tp] of 
the dynamics governed by H terminates the pulse. This 
motivates the two following ansatzc for the two unitary 
operators in Eq. [S] 

Up{Tp,Ts) = e-<^'^~^^'>" e-'^^^^^s U2{Tp,Ts) (7a) 

Up (t„ 0) = f/i (t„ 0) « fo' -Wdt ^-^r^H^ (7,^) 

which define the operators Ui{ts,0) and U2{Tp,Ts). 
Both, Ui and U2 must be seen as the corrections which 
are necessary because the left hand sides in Eqs. [7] do 
not equal the right hand sides without them since the 
pulse and the bath Hamiltonians do no commute. Note 
that no time-ordering is necessary for the exponential of 
H because it is constant. For the exponential of Hq no 
time-ordering is necessary despite the time dependence 
because the commutator [Ho{t), Ho{t')] — for any two 
times t and t' . 

If the corrections are small, e.g. U2{tp,Ts)Ui{ts,0) — 
1 -|- 0{Tp) up to linear order in Tp, this implies 



UpiTp,0) 



-i{Tp 



)H g-i<Ty Jo" vit)dt^- 



iTsH 



(8) 



which corresponds to the ideal, instantaneous pulse 
sketched in Fig. [T}) at Ts while the dynamics before and 
after Tg includes both the bath dynamics and the cou- 
pling between qubit and bath. 

The unitary corrections Ui and U2 are determined from 
the Schrodinger equation. We start analyzing U2 (I7ap . 
Let r be a time instant in the interval ]Ts,Tp] such that 



Up{T,Ts) = e — ■■- e 
where we use the difference 
Ax 



_ -iArH -iT Ho(t)dt 



U2iT,T,), (9) 



(10) 



generally for any variable x; here in particular for x = t. 
The Schrodinger equation of Up (t, r^) is the usual one 

idrUp (r, Ts) = (HoiT) + H) Up (r, r,) . (11) 

Inserting Eq. ^ into (fTT|) one obtains 

l^g-.ArH ^ e-'^^^^HoiT)} e-'^^^. U2 (r, r,) 

= (i/o(r) + H)e-''^^" e~'^r^ U2 (t,t,) . (12) 

It is easy to see that the terms with H cancel so that the 
differential equation for U2 becomes 



idrU2 (r, Ts) = ( Ho(t) - Ho{t) ) U2 (r, r,) , (13a) 
where 



^^o(r) 



— ^^/x" Ho{t)dt JAtH 



HoiT)e 



Formal integration leads to 



U2iT,Ts)^T 



exp —I 



-iAtH ^-i Hoit)dt 

(13b) 
(14a) 



F{t)dt 



where 



F{t) ■.= Ho{t)~Ho{t). 



(14b) 



From its defintion it is obvious that F{t) vanishes iden- 
tically if Hq and H commute. This is the case if there is 
no coupling between qubit and bath. Closer inspection 
of Eqs. (|13bll4b|) shows that F{t) = 0{tXv{t)). More 
generally, F(t) can be expanded in a series of H, i.e. in 
the parameters tX and tu which amounts up to an ex- 
pansion in Tp. Thus our approach provides the intended 
expansion in the shortness of the pulse. 
With the definitions 



ay (At) 



iAtH „ -lAtH 

e ayC 



(15) 



and ^ the time-dependent operator F{t) can be written 
in the compact form 



Fit) = vit) 



TyLp{t) 



(Tj,(Ai) e 



~iay^{t) 



(16) 
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Next, we treat Ui by the analogous procedure solving 
the Schrodinger equation for (ffb)) . Let r be a generic 
instant in the interval [0, rs[, the ansatz for this interval 
is 



Up (t„ r) = Ui (t„ r) e''^(^)'^« e^^^^. 



(17) 



To obtain the Schrodinger equation in the second time 
argument of Ui{ts,t) it is convenient to consider the 
hermitian conjugate of (llip yielding 

- idrUp (t„ r) = Up (r„ r) [FqW + H] (18) 

where the property U^ (^1,^2) = Up (t2,ii) has been ex- 
ploited. The procedure is the same as for U2 with the 
only difference that now the exponentials occur on the 
right side of Ui 



idrUl{Ts,T) = Ul{Ts,T)F{T) 



(19) 



with the same F{t) as defined in (I14bp . Formal integra- 
tion yields 



C/i(T„r) = T 



exp —I 



F{t)dt 



(20) 



Finally, we combine both corrections Ui and U2 

Up{Tp, 0) = Up {Tp, Ts) Up (Ts, 0) 



where 



Uf{tp,0) C/2(Tp,r,)[/i(T„0) 



T e 



}T{. 



(22) 



Thus the total correction Up is given by the time-ordered 
product of F(t). Thereby, we have derived the general 
expression for the difference between the unitary action 
of the real pulse as sketched in Fig. [T^ and the idealized, 
instantaneous pulse in Fig. [TJ). 



B. Expansion to Second Order in TpH 

We want to find the conditions under which Upijp^Q) 
can be approximated by the identity operator. Because 
F{t) is 0{X) for A — > the expansion in F is an ap- 
propriate first step. A convenient way to obtain the 71-th 
order of an expansion in A and loi, is to compute up to the 
71-th order of the Magnus expansion from average Hamil- 
tonian theory [§, [2^ . Then the resulting expressions are 
expanded to n-th order in A and w;,. The advantage of 
this approach over a direct expansion is that the expan- 
sion is done in the argument of the exponential. 



The Magnus expansion reads 

Uf{tp,Q) = T 
= exp 



exp [ -1 / F{t)dtj 



([-irp(F(i) + + + ...)], (23) 



where F^^^ — J^'' F{t)dt is the average value of F. 

The next-leading term F^^^ comprises the commutators 
of F with itself at different times 



Jo 



dh[F{t{),F{t2)]. 



(24) 



Thus from the Magnus expansion we know at least the 
two leading orders in powers of F{t). Next, we expand 
F{t) itself which is still a complicated quantity, see Eq. 
(Uni). To obtain the two leading orders in H we expand 
F{t) in powers of H corresponding to an expansion in Tp. 
We use the identity [22|] 



a,(Ai) = c7, + ^-(A0"[[i?,a,]]„ 



(25) 



where Ai is used as defined in pUj) . The notation stands 
for 

[[H,Gy]]i ^ [H,Gy\ ^ -2ia,XA 

[[H, <jy]]2 = [H, [H, cTy\] = -2ic7,X[Hb, A] + AayX^ 



(26) 



[[i/,a,]]3 = [i/, [H, [i/,a,]]] = ..., 
to obtain 

(Ty(At) - CJy = 2(Tj;AtAA "f 

{i(j^X[Hb, A] - 2ayX^A^) {Atf + 0{At^). (27) 

Inserting the above expansion in the equation (|16p of F(t) 
and using the elementary relation 



^i'^yvW g-i<T,v(t) = cos{^p)a^, + sin(^)(T, 



(28) 



where ip{t) = 2ip{t) as well as the Magnus expansion 
([23]) allows us to compute the leading orders of Up in 
exponential representation 



Uf{tp, 0) = exp[-i(?7^^^ + 77^^^ + 

with 

= iVll'^x + Vi2crz)XA 



(29) 



(30a) 



= i {mi(^x + ri22(^z) X[Hb, A] + ri23ayX^A\ (30b) 



Note that [Hh,A] = 0{uji,) so that 77^^^ truly represents 
the desired second order in the shortness of the pulse. 
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The coefficients rjij are given by the integrals 





l-Tp 

= 2 / Atv{t)cos{'4!{t))dt 
Jo 


(31a) 


Vl2 


= 2 Atv{t) sm{il;{t))dt 
Jo 


(31b) 


V21 


= / Afv{t)cos{ip{t))dt 
Jo 


(31c) 


V22 


= / At'^v{t)sm{ilj{t))dt 
Jo 


(31d) 


V23 


= -2 / {Atfv{t)dt 
Jo 


(31e) 



+4/ dti I dt2AtiAt2ViV2sm{'ipi - 1P2), 
Jo Jo 

where we used ipi for V-'(ij) a-nd Vi for v{ti) for i G {1,2} 
in the last line. The differences are used as defined in 
(fTU)l . The last line of piep results from the commutator 
([M)) yielding 

[F{ti),F{t2)] = AtiAt2'UlW2 ([CTz,Cr:r] sin V'l COSV'2 

+ [o'x, Cz] costAi sin 1/^2) A^A^ 
= 2iAti At2WiW20'j/ X 

(sin ipi cos ^02 — cos ipi sin ■02 ) 
= 2iAtiAt2ViV2cry sin (-01 - 1P2) A^A^ (32) 

This concludes the expansion up to second order in TpH, 
i.e., in the shortness of the pulse. 

IV. SHAPING THE PULSES 
A. General Discussion 

1. Linear Order 

The general idea is that one can shape the pulses such 
that the leading deviations in Up from unity vanish. 
First, we focus on the linear order which requires that 
rj^^^ vanishes. Hence we have two conditions given by 
rjii — and ryi2 = 0. The relation v{t) = dtilj{t)/2 can be 
exploited to integrate the corresponding integrals (|31ap 
and (|31bp by parts yielding 



= {tjj — Tg) sin ijj{Tp) + Ts sin ^{0) — sint/jdi (33a) 



= (rp — Ts) cos V-'(Tp) + Ts cos V'(O) — / cos-^di. (33b) 

Jo 

Any pulse which fulfills these two conditions will show 
only quadratic deviations from the idealized instanta- 
neous pulse. Note that Tj is not a priori fixed and can be 
considered a free parameter which can be tuned to fulfill 
the above conditions. Hence one additional free param- 
eter is sufficient to obtain a pulse which is ideal up to 
linear order. Explicit solutions will be discussed below. 



A last general property worth mentioning is the sym- 
metry of Eqs. (j33p under the transformation v{t) 
v{Tp—t) which implies — > Tp—Ts and ip{t) — > —'ip{Tp—t). 
Thus, ifv{t) fulfills Eqs. ^ for r^, then v{Tp-t) fulfills 
them for r„ — r, also. 



2. Quadratic Order 

The requirement 77^^^ in quadratic order adds another 
three integrals for 7721, ??22, and 7723. It is obvious that 
integration by parts can also be used to make v{t) disap- 
pear in the expressions for 7721 , 7722 , and 7723 . The resulting 
expressions are 



'721 = 



?722 



^^P—Iflsin^{rp)~^sinm 
At sin ij{t)dt 

(ZiL_Z£)!cosV(rp) + ^cosV(0) 



At cos ip{t)dt 



(34a) 



(34b) 



?723 = {Tp - Ts)TsSin0 - Ts / sm('0(i) - 7/;(0))rft 

-{Tp-Ts) / sin{-ip (Tp) - ^(t))dt 
Jo 

+ ^ J J sin(V'i - '02)sgn(ti - i2)rfiirfi2, (34c) 

where we used 6 = ipirp) — V'(O). 

In view of the five expressions to be set to zero, one 
needs at least five free parameters including Ts, in order 
to obtain pulses which show only cubic deviations from 
the idealized instantaneous pulse. One way to approach 
the solution of the above discussed conditions is to con- 
sider symmetric pulses aiming at Ts = Tp/2. Then ?/'(i) 
is an odd function about and the coefficients 7711 and 
7722 vanish by antisymmetry. Then only three conditions 
remain to be solved. But the shape of the pulse is already 
fixed to some extent by its symmetry. 

A closer inspection reveals that tt pulses cannot be 
corrected in second order. This statement is proven rig- 
orously by the following consideration. We combine the 
two real Eqs. p4ap and (|34b[) to the complex equation 

= r>''''(°) - (Tp - Tsfe'"^^^"^ + 2 / " Ate'^f'^'Ut, (35) 
which leads to 

^2 + „ ^^)2 ^ _2 r Ate:(^^'')-^'-°^^dt 



< 2 \At\dt 

Jo 

< r2 + (rp-T,)2. 



(36) 
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Figure 2: Composite tt and 7r/2 pulses with piecewise con- 
stant amplitude in the symmetric case. In the n pulse, the 
sign changes occur at ri = Tp/7 (6rp/7, respectively) and the 
amplitude is flmax ~ 7n /6[1/Tp]. In the tt/2 pulse, the sign 
changes occur at ri = 0.13155rp (r^ — ri, respectively) and 
the amplitude is amax ~ 1.65765[l/rp]. 

The equality holds if and only if the exponential factor 
in the first line of ((36)) changes its value from — 1 to 1 
abruptly at At = 0. This implies an instantaneous, ideal 
pulse with infinite amplitude which is not experimentally 
realizable. Hence a tt pulse cannot be corrected in second 
order. This no-go result is very remarkable in its general 
validity. 

For tt/2 pulses we found solutions which make 7721 and 

7722 vanish. In spite of our intensive search we have not 
succeeded in finding solutions which additionally make 

7723 vanish. Hence we conjecture that no such solution 
exists. But to our knowledge, there is no mathematical 
proof stating the impossibility to correct tt/2 pulses in 
second order. 

In this context, the reader may wonder whether one 
cannot combine two tt/2 pulses to obtain a tt pulse. But 
this is indeed not possible within the present framework. 
According to pT|) there is always some decoherent time 
evolution before and after the ideal pulse. Hence no di- 
rect combination of pulses with small angles to pulses 
with larger angles is possible. But pulse sequences with 
finite time intervals between the pulses as in dynamic 
decoupling are well possible [J, i, S, i, i, i, S H • 



3. Possibility of Ts — Tp ? 

We return to the corrections in linear order p3p . In- 
spired from discussions with experimentalists in NMR, 
we pose the question whether it is possible to shape a 
pulse such that it corresponds to an ideal, instantaneous 
pulse at the end of the real pulse, i.e. Tg = Tp. If such a 
pulse existed it could be used to measure the response of 



Figure 3: Composite tt and tv/2 pulses with piecewise constant 
amplitude made from two elementary pulses in the asymmet- 
ric case. For the tt pulse, the sign change occurs at ri — 3rp/4, 
the amplitude is flmax ~ 7r[l/rp] and the instant of the equiv- 
alent ideal pulse is Ts = (rp/2)(l — I/tt). For the 7r/2 pulse, 
ri = 0.78220rp, the amplitude IS dmax — 

1.39156[l/rp] and 
the instant of the equivalent ideal pulse is Ts — 0.23128rp. 

systems at delay times much shorter than the duration 
Tp of the initial pulse. Clearly, this would represent a 
promising experimental tool. 

Unfortunately, it is impossible to construct such a 
pulse. Again, we can prove this rigorously by inserting 
Ts = Tp into the equations p3p which we combine to one 
single complex equation using exp(z'0) = cos -\- i smip 

j-Tp 

= rpexp(ii/>(G)) - / exp(i?/>(t))dt. (37) 
"'0 

We find again a rigorous estimate 

Tp = ( \^v{i[^{t) - i'{0)]}dt 

Jo 

< I Idt = Tp (38) 
Jo 

which is fulfilled if and only if 7/;(t) is piecewise constant 
with jumps of 2tt. But such a ipit) represents via dttl}{t) = 
2v{t) instantaneous pulses with infinite amplitudes. Thus 
there cannot be a real pulse with bounded amplitude 
which is close to an ideal pulse at the end of its duration 
in linear order of the bath parameters. 

One can only approximate the desired situation by 
aiming at a close to Tp. To this end, asymmetric pulses 
will be considered below. 

B. Solutions 

In view of the above general discussion we focus here 
only on the correction in linear order. First, we consider 
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Figure 4: Composite tt and tt/2 pulses with piecewise con- 
stant amplitude made from three elementary pulses in the 
asymmetric case. For the n pulse, the sign changes oc- 
cur at ri = Tp/13 and at T2 = 6rp/13, the amplitude is 
tSmax ~ 137r/6[l/rp] and the instant of the equivalent ideal 
pulse is Ts — Tp/2. For the it/2 pulse, the sign changes occur 
at n = 0.44834rp and at T2 = 0.80813rp, the amplitude is 
amax = 2.80074[l/rp] and the instant of the equivalent ideal 
pulse is Ts = 0.219Tp. 

so-called composite pulses which consist of piecewise con- 
stant pulses of maximally positive or negative amplitude 
iimax- Second, continuously shaped pulses made from 
sines and cosines are investigated. In each class, we look 
at symmetric pulses with = Tp/2 first and then at 
asymmetric pulses. 

1. Composite Pulses 

a. Symmetric Pulses For symmetric pulses with 
Ts = Tp/2 there is only one equation to be solved, namely 
P3b[) . because p3ap is fulfilled by antisymmetry. We 
need only one free parameter and take the instant ri, at 
which the pulse changes first, to be this free parameter. 
Another sign change occurs by symmetry at t = Tp — ti. 
Fig. [5] depicts the solutions for a tt pulse and for a 7r/2 
pulse. The parameters are given in the caption. 

It is noteworthy that the symmetric tt pulse found is 
exactly the one that Cummins et al. named SCORPSE 
[23I [24 | . Do our results reproduce theirs generally? The 
answer is "no" because the goals are different. While 
Cummins et al. aim at finding a pulse shape such that 

Up{Tp,0)^Ue (39) 

our aim is to approximate as closely as possible 

Up{Tp, 0) « exp{-i{Tp - Ts)H)Ug exp{-iTsH), (40) 

where Ug stands for the ideal, 9 pulse evolution operator. 
In other words, the authors of the Refs. [H, [2J aim at 



Figure 5: Continuous tt and 7r/2 pulses in the symmetric case. 
The paramter a as defined in the ansatz (|42p takes the value 
— 3.73002/rp for the tt pulse and the value — 1.78259/rp for 
the 7r/2 pulse. The horizontal lines indicate the amplitude 
for a totally constant, uncorrected tt pulse (solid line) and 
7r/2 pulse (dashed-dotted line), respectively. 



making the imperfections vanish completely while our 
aim is more modest: we aim at separating the bath and 
the pulse evolution. We are convinced that our aim is 
more realistic, in particular for general dynamic baths 
without static components so that no averaging to zero 
is possible. The compensation of the coupling between 
qubit and bath for longer storage times shall be done by 
dynamic decoupling [J [3, H, 0, IE B H Hi on the basis of 
the optimized pulses discussed here. 

But for static couplings H (x and a tt pulse with 
Ts — Tp/2 both aims p9l40p happen to coincide 

= U^- (41) 

Note that the argument does not require that the pulse 
is symmetric but only that the splitting time Ts is in the 
middle of the pulse so that the first half of the evolu- 
tion under H compensates the second half. Indeed, the 
CORPSE pulse [Hill fulfills our conditions §^ as weU, 
see Fig. g] 

b. Asymmetric Pulses On considering asymmetric 
pulses there is no reason to restrict oneself to = Tp/2. 
Hence there is another free parameter to tune. In this 
way, it is possible to find composite pulses which consist 
only of two constant regions. Fig. [3] depicts the simplest 
solutions, but there can be many more. For the tt pulse, 
for instance, we have found an infinite set of solutions 
with the time ri = (2n 4- l)Tp/{An) at which the jump 
occurs and the amplitude Omax = nn/Tp where n is a 
positive integer. The corresponding splitting time reads 
Ts/Tp = 1/2 + (— l)"/(2n7r). The solution depicted cor- 
responds to n = 1. 



The pulses depicted in Fig. [3] are the most asymmetric 
ones for given maximum ampUtude. This is reflected by 
the values for r, which are not close to Tp/2. If one is 
looking for a pulse shape which allows the measurement 
of a signal as soon as possible after an ideal pulse the 
inverted pulses v{Tp — t) with Tg ~ 0.77 for the 7r/2 pulse 
and with « 0.66 for the tt pulse should be used. 

Beyond the composite pulses consisting of two con- 
stant regions we consider also the more commonly dis- 
cussed composite pulses consisting of three constant re- 
gions. Fig. m shows generic results. The tt pulse depicted 
is the CORPSE pulse proposed previously [H, H^l which 
happens to fulfill also our conditions as explained above. 
The solution is arbitrary in the sense that other values 
for Ts are also possible. 

The same is true for the 7r/2 pulse for which we have 
chosen at will. Many other values are also possible. 
For Ts = Tp/2, however, we have not found any asym- 
metric solution. 



2. Continuous Pulses 

The above composite pulses are advantageous because 
they are fairly simple to generate and because they use 
the maximally possible amplitudes most efiiciently. But 
they are not optimum in view of their frequency selectiv- 
ity due to the jumps. It is well-known that the faithful 
representation of jumps requires particularly broad fre- 
quency bands. Hence it is important to demonstrate that 
our conditions (|33p can be also fulfilled by continuous 
pulses. We choose ansatze inspired from Fourier series 
which ensure even continuity of the derivative v' {t). 

a. Symmetric Pulses For symmetric pulses our 
ansatz reads 

v{t) = 9/2 + a cos(27rt/Tp) - (a + 6'/2) cos(47ri/Tp). (42) 

The coefficients are chosen such that u(0) = v{Tp) = 
and w'(0) = v'{Tp) = 0. Of course, higher cosine terms 
could be added. But the constraints in linear order ([55]) 
can be fulfilled already by the above ansatz. Fig. O dis- 
plays the resulting pulse shapes. Note that the instant of 
the equivalent ideal pulse is = Tp/2. Comparing Fig. 
[5] and Fig. [5] the above statement is illustrated that the 
composite pulses require only smaller amplitudes. But 
the continuous pulses in Fig. [3] are much smoother and, 
hence, do not spread so much in frequency space. It de- 
pends on the actual constraints in experiment which kind 
of pulse is the most advantageous. 

6. Asymmetric Pulses For completeness, we also 
consider an ansatz allowing for asymmetric continuous 
pulses 

v{t) = {0/2){l-cos{2-Kt/Tp)) 

+afiin{2TTt/Tp) - (a/2) sin(47rt/rp). (43) 

The resulting solutions of the conditions ([55)1 are plot- 
ted in Fig. [51 The values of the parameters are given in 
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Figure 6: Continuous tt and 7r/2 pulses in the asymmetric 
case. For the tt pulse, the parameter a as defined in the 
ansatz (|43p takes the value 3.39308/rp and the position Ts of 
the equivalent ideal pulse takes the value 0.33806rp. For the 
n/2 pulse, the parameter a is 1.54714/rp and the position Ts 
of the equivalent ideal pulse is 0.28093rp. The horizontal lines 
indicate the amplitude for a totally constant, uncorrected tt 
pulse (solid line) and 7r/2 pulse (dashed-dotted line), respec- 
tively. 

the caption. Naturally, the instant of the equivalent in- 
stantaneous pulse is no longer in the middle of the pulse 

Ts ^ Tp. 

V. CONCLUSIONS 

In this paper, we have investigated under which con- 
ditions it is possible to shape the pulses implementing 
single qubit gates such that they correspond to ideal, in- 
stantaneous pulses (cf. Fig.[T]). To this end, we computed 
the corrections in powers of H, the Hamiltonian contain- 
ing the coupling to the bath and the dynamics of the 
bath, to the equivalence 

Up{Tp, 0) = e-'^^'--'-'^"Uee-"'" + 0{H'^). (44) 

The expansion in powers of H is an expansion in the 
shortness of the total pulse. No assumption about the 
nature of the bath is made so that our results are very 
generally applicable. 

We have derived the explicit expression for the cor- 
rections in linear order in H as well as those quadratic 
in H . In the evaluation of the expressions found we fo- 
cussed on the tt pulse and the 7r/2 pulse. They are by 
far the most important ones for all kinds of applications. 
The TT pulse is used to flip spins or qubits. An important 
application arc spin echo experiments and, more gener- 
ally, dynamic decoupling of two-level systems from their 
environment [ll, |3|, U, i, i, S i] ■ The 7r/2 pulse gener- 
ates in NMR experiments the precessing magnetic field 
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and is used in many intricate pulse sequences to suppress 
the effect of unwanted interactions Q . In the context of 
quantum information, the 7r/2 pulse realizes the impor- 
tant Hadamard gate. 

We found and showed a multitude of pulse shapes 
which are correct in linear order, i.e. 

Up{Tp, 0) = e-''("--^=)^[/ee-''"=^ + OiH^). (45) 

The pulses can be chosen piecewise constant which re- 
quires the lowest amplitudes. They imply, however, 
jumps which deteriorate the frequency selectivity [13] ■ 
Continuous and continuously differentiable pulses can 
also be realized and they have a much better frequency 
selectivity. Their drawback is that the maximum ampli- 
tude required is larger than for the piecewise constant 
pulses. 

We noted that our objective (1451) ha ppens to coincide 
with the one of earlier work [l^, I23ll24| for a tt pulse with 
Ts = 1/2 applied to a completely static bath, i.e., a bath 
without internal dynamics. One reason for this coinci- 
dence is that in leading order the internal bath dynamics 
does not play a role. The commutators between different 
bath operators occur only in second order. Hence the cor- 
responding CORPSE and SCORPSE pulse fulfih also the 
conditions derived in the present paper. In general, how- 
ever, the coincidence does not hold, i.e. if Ts ^ Tp/2 or if 
the bath is not static in the above sense or if pulses with 
angles 9 different from tt are considered. Moreover, we 
emphasize that for dynamic baths without static com- 
ponents, i.e., baths with internal dynamics, the aim to 
make the coupling between qubit and bath vanish alto- 
gether, i.e. Up{Tp,0) « Ug, cannot be fulfilled for general 
angle 9 ^ tt. 

Investigating the quadratic order, we could prove that 
TT pulses cannot be corrected in this order. For 7r/2 
pulses we have not found any solution in spite of inten- 
sive search. Hence we are led to conjecture that no such 
solution exists. It would be interesting to find a math- 
ematical foundation for this conjecture. Even in linear 
order, it is not possible to shape the pulse in such a way 
that it corresponds to an ideal, instantaneous pulse at 
the very end of the real pulse. This would have been a 
nice property for NMR measurements since it would have 
permitted to measure directly after a given pulse with- 



out further delay. One may, however, use asymmetric 
pulses which correspond to ideal pulses at about 77% of 
the total pulse length. 

The message from our findings is that pulses can be 
shaped such that they approximate ideal, instantaneous 
pulses. The required pulse shapes have to fulfill rather 
simple analytic integrals so that everyone can fine-tune 
his pulse shape easily himself. The optimized tt pulses 
are an excellent starting point for optimized dynamic de- 
coupling schemes [8|. 

Mottonen et al. [101 have posed themselves a simi- 
lar question for tt pulses as we have done. The main 
differences are that they investigated classical noise nu- 
merically while we tackled a fully quantum mechanical 
model by analytical means. Another difference is that 
they aimed at Up{Tp,0) ~ C/^ whereas our goal was not 
to make the bath vanish, but to disentangle it from the 
actual pulse, i.e., the coupling between qubit and bath 
is still present before and after the idealized pulse. Of 
course, in a model of classical noise there are no bath op- 
erators which do not commute so the mere disentangle- 
ment is trivially given in the model studied by Mottonen 
et al.. We found that Up{Tp,0) « U-^ and Up{Tp,0) = 
^~i{Tj,/2)Hjj^^-i{Tp/2)H coincide for baths without inter- 
nal dynamics. For these special cases, our pulse shapes 
agree with those previously found [lo, [131 ■ For in- 
stance, the SCORPSE pulse appears to be a very good 
choice. 

Further work should simulate the proposed pulse 
shapes in specific models in order to understand better 
how important the neglected corrections in second order 
really are. It is to be expected that for small values of A 
(generic coupling to the bath) and LUf, (generic frequency 
of the bath) the corrected pulses are superior to the sim- 
ple, uncorrected ones. But if the characteristic times of 
the bath are too short, presumably the simple pulses will 
be better. 
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